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We discuss the necessity of non-perturbative renormalization in QCD and HQET and explain the general 
strategy for solving this problem. A few selected topics are discussed in some detail, namely the importance of 
off-shell improvement in the MOM-scheme on the lattice, recent progress in the implementation of finite volume 
schemes and then particular emphasis is put on the recent idea to carry out a non-perturbative renormalization 
of the Heavy Quark Effective Theory (HQET). 



1. INTRODUCTION 

Non-perturbative renormalization is an impor- 
tant mile-stone on the path from lattice QCD 
computations to precise predictions for particle 
physics phenomenology. Why is this so? One 
might think that only the QCD coupling and 
quark masses need to be renormalized. However, 
very important phenomenological applications of 
lattice QCD concern the physics of heavy quarks 
and weak decays. As illustrated in Fig. |^, these 
involve energy scales (masses), which are far too 
high for the lattices that can be simulated on 
present and (near future) computers. Effective 
theories have to be used, summarizing the effects 
of the heavy fields, which can't be treated dy- 
namically, in local composite operators. They 
are a way to implement expansions in terms of 
1/Mx, X = W, t,.... The requirement that the 
corrections to the n-th order are really (1/Afj:)"+^ 
and terms such as as/(Afa;)" do not appear fixes 
the renormalization of the operators in the effec- 
tive Lagrangian. We denote the so-defined renor- 
malization scheme by "match" . 

To give a specific example, the mixing am- 
plitude between a K-meson and a K is to low- 
est order in the weak interaction given by a 
QCD matrix element of a 4-fermion operator 
whose (re)normalization is completely fixed by 
the matching to the standard model amplitude. 

Of course, already before considering weak de- 
cays, the renormalization of the QCD coupling 
and quark masses is important to obtain the ba- 
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sic renormalization group invariants (RGI) of the 
theory: the Lambda parameter, Aj^, and the 
RGI quark masses, Mi,i = u, d, .... Their def- 
inition and their special role in parameterizing 
QCD is briefly but thoroughly explained in Q. 

Why should we do renormalization non- 
perturbatively? There are good reasons: 

1 . Concerning , Mi , phenomenological deter- 
minations claim very good precision but are fre- 
quently based on difficult to test assumptions, 
most notably the applicability of perturbation 
theory in various kinematical regions. In my opin- 
ion, the role of lattice QCD is complementary: to 
use as few assumptions as possible. This means 
to perform renormalization non-perturbatively. 

2. Some examples demonstrate that perturbative 
renormalization is not sufficiently precise in prac- 
tice. E.g. the JLQCD collaboration computed 
i?K with two different perturbatively renormal- 
ized operators which should be equivalent in the 
continuum limit. Performing the continuum ex- 
trapolations, they could only find agreement by 
fitting o? corrections in addition to the a} lat- 
tice artifacts A more recent example is the 
static-light axial current. As shown in Fig. ^ per- 
turbative and non-perturbative Z-factors differ by 
more than the unavoidable 0(a) term. 



2. GENERAL STRATEGY 

We now discuss the implementation of scale- 
dependent renormalizations, leaving aside simpler 
cases such as the renormalization of the axial cur- 
rent 0], where the (scale independent) renormal- 
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Figure 1. Top: The standard model summarizes physics at many different scales, while lattice QCD simulated 
by Monte Carlo (MC) can correctly cover only a small range between the inverse lattice spacing , and the box 
size L. Bottom: The strategy to connect the bare matrix elements to the ones in the matching scheme. 
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Figure 2. Non-perturbative (NP) renormalization 
constant of the static-light axial current |^ and its 
estimate from (tadpole-improved Q) perturbation 
theory. The label "improved" refers to the 0(a)- 
improved theory. 



ization may be determined from a symmetry. Fo- 
cusing on multiplicative renormalization, the bare 
operators, Obarc, and the renormalized ones at 
scale ^, Ointci (M), are related by 



CintoiV) = 2'intcr(go, aM)Cba 



(1) 



with a renormalization factor Z (or a matrix) de- 
pending on the bare coupling, go, and the com- 
bination a/Lt. In an intermediate renormalization 
scheme, Z is fixed by requiring 



'\Ointm■{^J■)\a) = (/3|C'bare|a)tr 



cclcvel 



(2) 



for convenient states |a), To have simple 

renormalization group equations (RGE) , it is im- 
portant that 1/3) are characterized through 
one scale, /i, only. In the two most frequently 



used intermediate schemes, this is realized as fol- 
lows. 

SF: In the Schrodinger functional schemes 0, 
eq. (H) is formulated through the gauge invari- 
ant path integral in a finite space-time volume 
T X L X L X L with T/L fixed. \a), |/3) are given 
in terms of some boundary states (propagated in 
Euclidean time) and one has ^ = 1/L. 
MOM [||: Here one uses infinite volume quark 
states in Landau gauge, with four-momenta p of 
the quarks satisfying = M^- 

Both are massless schemes, their renormaliza- 
tion constants are evaluated at zero quark masses. 
In addition, the renormalized operators are inde- 
pendent of the regularization.0 More precisely, 
(the limit a — > of) the matrix elements 



,{p) = (/lOintcr(M)N) 



(3) 



are independent of the regularization and the per- 
turbative coefficients 7^ in the RGE 



d$i, 



inter 



7(5^) = -[705^ + 71/ + .. 



(4) 
(5) 



can in principle be computed using dimensional 
or a lattice regularization. 

For phenomenological applications we are in- 
terested in the matrix elements 3> match (m) i'^ 
the matching scheme mentioned in Sect. |l|. To 
connect to it, it is usually]^ convenient to first 

^Thus also the name RI/MOM is used instead of MOM. 
^For an alternative see Sect. 3.3. 
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Figure 3. Typical relative error terms as they show up in infinite volume schemes (left) and finite volume schemes 
(right): a'^-efFects as well as non-perturbative (dotted) and perturbative (dashed) terms. 



compute the scheme independent renormalization 
group invariant (RGI) 



$RGi = hm $i„tcr(/^) [26o5^(m)] 

/J.— >-oo 



^intcr(M)(26o.9(M) ) 

X exp 



(6) 
(7) 



and then connect $rgi to $match i^J■) through the 
analogue of eqs.(^,^) for the matching scheme. 
This last step is mostly done by perturbation the- 
ory and - as discussed by P. Hasenfratz at last 
year's conference |9) - its reliability needs to be 
investigated. For the purpose of this report, we 
note that the RGI matrix elements eq. (^ are the 
fundamental quantities of QCD and should be de- 
termined from lattice QCD with good precision. 

The strategy is sketched in the bottom part of 
Fig. 0. For multiplicative renormalization, each 
arrow in this graph is realized by a factor. 



match 
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Z inter (go, ^J■ma) X $barc(5o)- 



(8) 



The first two factors are independent of the lattice 
regularization (e.g. choice of action). 

2.1. Challenges 

The practical implementation of this pro- 
gramme represents a challenge beyond the one 
present in spectrum calculations, since the limit 
At —> 00 in eq. (^ has to be controlled. 

With MOM (or other infinite volume schemes) 
as an intermediate scheme, one computes 



*inter(A<) with Zi^terigo, a^) from eqs.(|]J|) for ^ = 
0(2 GeV) and uses the perturbative RG eq. (||) 
inserted into eq. (^ to continue to /i ^ oo. Sys- 
tematic errors of the order shown on the l.h.s. of 
Fig. § have to be expected. It is evidently not 
easy to disentangle the various error sources in 
the window of /i and a available on large volume 
lattices. 

To cleanly separate the lattice artifacts (O(a^) 
or 0(a)) from the physical /z-dependence, renor- 
malization conditions may also be posed in finite 
volume (e.g. SF-scheme) and then a recursion 
^ ^ 2^ allows to reach large enough /i to make 
sure that the terms illustrated on the r.h.s. of 
Fig. H are negligible/controlled pO| , pd| . The key 
point is that before comparing the /i-dependence 
to the perturbative one, the a-effects can be con- 
trolled in each step of the recursion. Examples 
are given in Sect. |3.2|. 



3. SELECTED TOPICS 

We now discuss some selected topics of rele- 
vance. Recent reviews, with partly different em- 
phasis, are JT2| . 

3.1. MOM-scheme 

The MOM-scheme [|| is very popular. It is 
rather easy to implement and does not require 
simulations to be done specifically for the renor- 
malization. A number of new numerical results 
have been presented at this conference ||l^,|lj,|l5j . 
Here we concentrate on two aspects which have 
not been discussed much. 
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3.1.1. O(a)-improvement 

Renormalization conditions in the MOM- 
scheme are based on off-shell quark Greens- 
functions in Landau gauge, which are not 0(a)- 
improved by the standard on-shell improvement 
programme |]l6| , pT| . The most basic ingredient is 
the momentum space quark propagator 

~S{p) = J2e~'P^^Pix)m) (9) 

X 

= [^7^p^I]l + . (10) 

For sufhciently large , it is computable in per- 
turbation theory and S2 is expected to go to a 
constant (asymptotic freedom!). Evaluating E2 
with the non-perturbatively on-shell improved ac- 
tion (e.g. at a = 0.07 fm), one finds in- 
stead a strong rise with p"^ (• in Fig. |^). 
The SPQcdR Collaboration verified that this lat- 
tice artifact is reduced for decreasing a jlj]. We 
learn that off-shell 0(a) lattice artifacts are large 
and their impact on the numerical results in the 
MOM-scheme may be substantial. 

The possibilities of reducing them in a sys- 
tematic Symanzik improvement programme were 
studied in [po|j2l|] . A rough summary follows. 

Since one is interested in non gauge invariant 
Greens functions, the improvement terms have to 
satisfy only BRST invariance instead of gauge in- 
variance. In addition, contact terms have to be 
subtracted from the correlation functions in po- 
sition space. 

We restrict ourselves to the quark propa- 
gator and correlation functions of the form 
(V'(a;)V^(y)C'r(z)), with Or = i^TT''^P, T a ma- 
trix in Dirac-space and a traceless matrix r" in 
flavor space. These are the Green's functions 
needed for the MOM-renormalization of flavor 
non-singlet 2-fermion operators. Then off-shell 
0(a)-improvement may be achieved by pO[ | 

1) replacing the quark fleld ip^x) with 

i;R{x) = Z-^/^{l- ^b^amq)Mx) (H) 
V'i(a;) = [1 + ac'^D + acNGid^t1t,]tpix) , 

with (7o-dependent coefficients Zq, 6q, c^, cngi 
(and similarly for ijj). 

2) adding a term acYtpTDtp to Or in addition to 
the on-shell improvement terms of |16[ . 



Above, D denotes the full lattice Dirac opera- 
tor including the mass term. An important point 
to note is that this is strictly speaking not an 
improvement of the fields: e.g. an A''-point func- 
tion of Or{x) will be 0(a)-improved when consid- 
ered on-shell (all operators separated by a phys- 
ical distance) but not when considered off-shell 
(e.g. Fourier transformed). In the latter case, 
additional contact terms have to be subtracted 
in general. Nevertheless, the above is a possible 
form of writing the 0(a) counter- terms needed for 
the MOM-scheme renormalization of Or . 

Investigations of practical ways to determine 
the new coefficients Cq,CNGi have started [ p^|j2^ . 
Before we discuss the result, we turn to off-shell 
improvement in the context of formulations of 
QCD with exact chiral symmetry (see Q). Their 
massless Dirac operator D satisfies the Ginsparg- 
Wilson relation Q 

j5D-^ + D-^j5=a2Rj5, (12) 

with a local operator R and the fermionic action 
is S = a^ + moil ~ ^RD)]ij{x). 

Examples of local Dirac operators satisfying 
eq. ([1^ ) are known [p5|j2^ and in addition it has 
been shown that Domain wall (DW) fermions ||2^ 
are approximate realizations ||28| . 

A possibly confusing point is that the MOM- 
scheme is applied for these formulations, stating 
that there is 0(a)-improvement because of chiral 
symmetry, while Liischer's exact chiral symme- 
try is only valid on-shell (otherwise it would 
be in contradiction to the Nielsen-Ninomia theo- 
rem |3^). We briefly explain in which sense off- 
shell 0(a)-improvement is implied by eq. (|l2|); see 
PH for a similar discussion. For simplicity we as- 
sume that R is proportional to unity. 

Introduce quark fields and anti-quark fields 

M^) = - ^RDMx) , ^iix)=^{x), (13) 

and perform the same replacement in quark bi- 
linears, 

Or,i(a:) =^i(a;)rT>i(a;). (14) 
It follows immediately that the quark propagator 

Slip) = ^^"'"'(^1(^)^1(0)) (15) 



5 



2 - 



O 



staggered 
a = 0.12fm 

impr. Wilson 
a=0.07fm 

A DW, a = 0.1fm 

overlap 
° a=0.1frn 



10 



15 



20 (p r„)= 



1.5 



• impr. Wilson 
a = 0.07fm 



• . . * 



0.5 



I oog 



10 



15 



20 



(p 



Figure 4. The functions Ej, eq. (|lO[), renormalized such that they coincide at T,i{rQp^ = 5). Results are from 
(overlap), (DW), Q (improved Wilson) and |Q (staggered). Open circles use estimates of Cq,CNGi ||l9t , 
while filled ones are for c'„ — cngi ~ 0. Error bars are comparable to symbol sizes. 



satisfies 5^i(p)75 + 755i(p) = for mo = 0; the 
continuum relation is exact at finite lattice spac- 
ing. Furthermore, the Green's functions of quark 
bilinears, eq. (p^), satisfy the continuum chiral 
Ward identities at finite lattice spacing. This ex- 
cludes the necessity of additional 0(a) improve- 
ment terms since all possible terms [ pO| transform 
differently under chiral rotations and would lead 
to violations of the Ward identities. One con- 
cludes that the off-shell improved quark fields and 
quark bilinears are known without any free im- 
provement coefficients. 

Note that the action is not written in terms 
of the fields tpi, V'l since this would lead to a 
non-local Dirac operator (in agreement with the 
Nielsen-Ninomia theorem); only the fields in the 
correlation functions are improved. Concerning 
DW fermions in the limit of large "fifth dimen- 
sion" , the situation was clarified by Kikukawa and 
Noguchi. They derived a 4-d effective action with 
fields ip, ip and a local Dirac operator which sat- 
isfies eq. (12) with R — 2. The Shamir- Furman 
fermion fields ||2^, which are constructed from the 
5-d fields at the boundaries, are identified with ipi, 
■0J exactly related to ip, tp by eq. (p^. 

We now come back to numerical results. In 
Fig. ^ Si and E2 are displayed for various dis- 
cretizations. The results for DW fermions 
and overlap fermions show the expected lev- 



eling off of S2 with and are in complete agree- 
ment with each other within their (not so small) 
statistical errors. As discussed before, 0(a)- 
errors are expected to be absent^. The mu- 
tual agreement is an indication that 0(a^)-errors 
are small. More surprisingly, also the approx- 
imately off-shell 0(a)-improved Wilson fermion 
results [|l9) fit in very well! On the other hand, 
the staggered fermion E2 @1 appears to have sig- 
nificant discretization errors at a = 0.12 fm. 

What does this mean for NP renormalization 
in the MOM-scheme? Since the renormalization 
conditions, eq. (^), are imposed in such a way 
that Z = 1 + 0(^0) holds exactly at finite a, dis- 
cretization errors are suppressed by one order of 
(?o- Whether this is sufficient needs more inves- 
tigations. In any case, formulations with exact 
chiral symmetry [^2|j3^ are expected to be supe- 
rior and first numerical results for MOM-scheme 
Z's [^,^ appear to confirm this. 

At this conference, new results for MOM renor- 
malization constants in the on-shell improved the- 
ory were presented p^ , p^ . Not all of them are 
well compatible with the perturbative RG to the 
NL order. It remains to be clarified why this 
is so. Apart from a true non-perturbative /i- 
dependence, 0(a) errors or Goldstone pole con- 



* We assume that the extent of the fifth dimension in the 
DW computation was large enough. 
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tributions are possible explanations. 

3.1.2. Goldstone poles 

To discuss the latter possibility, we choose the 
simplest case, Os = ipl^T""^ ^^^d write schemati- 
cally 



R{m) = (p'lOsb) / (p'\0,\p) trcclcvcl ; 



(16) 



with (p'l {\p)) an (anti-) quark state. It is stan- 
dard to choose the forward kinematics p — p' and 
Zp^ in the MOM-scheme is to be obtained from 
R{m) in the limit where the quark mass m goes to 
zero. However, it is well known (see e.g. [pp^) 
that R{m) is singular in the chiral limit due to 
the coupling of the pion to the pseudo-scalar ver- 
tex and the fact that p — p' = {) in the chosen 
kinematics; one expects 



Rim) = A{p^) 



+ B{p') + 0{ml). (17) 



Based on this form it has become practice |13 TJ] 
to eliminate the first term by forming 



Bip') = 



ml{mi)R{mi) - m'^{m2)R{m2) 



+0{ml), 



(18) 



and then to determine Zp^ from the chiral limit of 
eq. (|l8|). From a theoretical point of view, eq. ( p7| ) 
means that the renormalization constant Zp (and 
others) do not exist in the massless MOM-scheme 
for p = p' . As alternatives, one could formal- 
ize the subtraction eq. ( p^ and include it into 
the definition of the renormalization scheme or 
one could choose a kinematics which truly corre- 
sponds to short distances in position space, i.e. 

(p')^j {p^p'Y ^ ^QCD o'" o^'S could avoid op- 
erators which couple to the pion and relate their 
renormalization to others by Ward identities |36|. 
The second alternative would certainly be a the- 
oretically sound solution. 

Concerning the numerical results alluded to 
above, the fact that one is dealing with the sub- 
traction of singular terms makes the chiral extrap- 
olation of the remainder numerically more diffi- 
cult. In fact, double poles have to be subtracted 
from the 4-fermion vertex functions f|l3|J37tl . 



3.2. Progress using finite volume schemes 

Finite volume schemes allow to compute the 
scale-dependence of Winter (/^) explicitly up to 
large = 0(100 GeV). One can then ver- 
ify the onset of perturbative running and use 
the perturbative (3 and 7 functions (in eq. (Q)) 
only a safe distance beyond that point. In this 
way one computes $RGi/*intor(Mm), cq. (||), non- 
perturbatively. Apart from the early studies 
of renormalized couplings JTop^ , only schemes 
based on the Schrodinger functional |^ have 
been used so far. The reasons are probably that 
in this framework 

• gauge invariance is explicit, 

• on-shell improvement is sufficient, 

• observables show good signal/noise ratios, in 
particular the important running coupling, 

• one can use a massless scheme, without 
extrapolations, 

• perturbation theory is relatively easy and 

• a-effects are typically quite small. 

Still, for new applications one should keep in 
mind that other finite volume schemes may be 
useful. 

The central objects, which describe the p,- 
dependence and are computable by MC-methods 
are the step scaling functions. They give the 
renormalized quantities at scale p = l/{sx L) 
as a function of those at scale = 1/L. (In this 
section we use just "R" to denote the intermediate 
scheme). Picking a complete set of observables, 
j, the step scaling functions F are 

<^>n,Ap/s)^F,{{^nA^')}) ' i,J = 0,...,Af .(19) 

A special role is played by the renormalized cou- 
pling 



$R,o = .9^(i), 



(20) 



which just like in eq. (^ is taken to parameterize 
the scale /i = 1/L. In fact, in a massless scheme, 
its step scaling function 

Fo=a{u) = g\sL)\ 



(21) 



needs no further argument. Another example of a 
step scaling function is the case of composite op- 
erators which mix under renormalization, where 
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eq. (|l9|) is realized as 



$R,,(/i/s) = ^a,,(u)$R,,(Ai) ,u = .g'(i) (22) 



^^,({$R,,.(^)})). 



additive renormalization see Sect. 3.3 



For an example of purely 
eq. 



3)- 

Lattice approximants of the step scaling func- 
tions are defined at finite resolution a/ L via 



a/ L) 



T.ij{u,a/L) = 



f{sL) 

at f{L) 



(23) 
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Figure 5. Step scaling function for Nf — 2 (circles) 
compared to the continuum a for A^'f = (squares). 



As indicated in the illustration for s = 2, the r.h.s. 
of eq. (|23| ) involves quantities on lattices L/a 
and sL/a but at the same bare parameters g^, mp 
which are fixed by (f' (L) = u and the vanishing of 
the PCAC-mass 111. Thus the lattice spacing a 
is kept fixed (and the quark mass zero) while the 
lattice size L changes by a factor s. In the next 
step one changes L/a at constant (downwards 
arrows in the graph), thus keeping L fixed in- 
stead. This process is iterated to connect L — Lq 
to L = Ln — s"Lo. To approach the continuum 
limit, all step scaling functions are computed for 
several values of a/L at fixed u and extrapolated 
to the continuum: (t{u) ~ lim^j/^^p S(m, a/L). 



/ 



Inverting the continuum step scaling functions 
allows to climb up in energy, /i ^ s/i — > s^^ .... 
Let us illustrate recent progress in implementa- 
tions of the approach. 

1. The basis for all subsequent applications is the 
running coupling. Following very closely |^,|ll], 
it has been studied also for N{ — 2. In addition 
to last year's results also T,{u,a/L = 1/8) 
is now available |Q. The lattice spacing depen- 
dence is very weak and the small effect of dy- 
namical fermions on the (SF-scheme) /^-function 



IS seen (Fig.|). Ahhough some of the simulations 
were difficult and a modification of the impor- 
tance sampling was needed to become convinced 
that the configuration space is sampled properly 
in all cases, one finally has a precision result for 
massless dynamical fermions! Compared to 41 
the systematic error due to the a/L — > extrap- 
olation is reduced, leading to 



-t'maxAjyjg 



0.68(7), where g2(L^^^) = 5. (24) 



2. Defining the quark mass through the renor- 
malized PCAC relation, its scale dependence is 
given by that of the renormalization constant, Zp, 
of the fiavor non-singlet pseudo-scalar density. In 
close analogy to Q, its step scaling function, 
Ep, has now been computed also for Nf = 2 and 
L/a = 6,8 [|ll. A continuum limit requires at 
least one larger value of L/a. 
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Figure 6. Running of quark masses for N{ = 2. For 
details see Q. 
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As a preliminary step, the running quark mass 
in units of the RGI-mass M was evaluated from 
Ep(u, 1/6) and Ep(w, 1/8) separately. The en- 
couraging result is shown in Fig. ||. Although 
the overall cost of the simulations grows very 
rapidly with L/a ||l[, using Hasenbusch's variant 
of HMC leads to an effort of only about a day 
on an APEmille crate (peak speed of 64 Gflop/s) 
for 1% precision in Z-p and for L/a — 16. Thus, 
the step L/a = 12 ^ 2iy/a = 24 is not expected 
to be a problem and the continuum limit is within 
reach. 

Assuming Aj^ = 0(250 MeV) one finds that 
the maximum box size, L,„ax, is large enough to 
make contact with the infinite volume physics and 
extract Aj^/F^ and Ms/F^. Still, performing a 
continuum limit in this step may be a challenge 
for the future. 

3. There are first results for the renormaliza- 
tion of parity-odd 4-fermion operators in the SF- 
scheme (quenched approximation) ||4^ . These 
are of considerable interest since their matrix el- 
ements give predictions for 

• parity-odd decays when their hadronic matrix 
elements are evaluated with the standard (im- 
proved) Wilson action 

• parity-even decays when the standard mass 
term is replaced by the "twisted" one see 
also @. We refer to this as tmQCD. 

The second option is particularly interesting since 
in this way Sk can be computed without any op- 
erator mixing.^ 

For the definition of the operators and their 
renormalization conditions, we refer to [Q. A 
complete basis has been considered, but so far 
the detailed analysis has only been done for the 
operator whose matrix element in tmQCD yields 
Bk- 

Although in the range a/L — 1/6 — 1/16 and 
for 13 different values of m a lattice spacing depen- 
dence of at most 5% has been observed in the step 
scaling functions, a continuum extrapolation has 
to be done with care; the operator is not 0(a)- 

^With the Wilson action and standard mass term, the rel- 
evant 4-fermion operator mixes with operators of wrong 
chirality. This has been a significant source of uncer- 
tainty ji^. 



improved and linear and quadratic a-effects may 
compete (see also the case of structure functions 
p8|| ). Fortunately there are a number of handles 
to control the continuum limit, a) Five differ- 
ent intermediate schemes have been implemented, 
which all should yield the same final RGI matrix 
element $r,gi- b) For the same scheme, different 
regularizations (e.g. Wilson and improved) have 
to extrapolate to the same continuum step scaling 
function [Q. These constraints should be suffi- 
cient to obtain accurate final numbers. 
4. The renormalization of the static-light ax- 
ial current in the quenched approximation has 
been finished, including nice continuum extrap- 
olations of the step scaling functions [Q. The 
overall Z-factor, mentioned already in the intro- 
duction, has then been used to estimate Fb^ = 
261(46) MeV -I- 0(l/mBj from the bare matrix 
elements of [^,^ in the unimproved theory. 

3.3. Heavy Quark Effective Theory 

As the b-quark mass is usually beyond the 
cutoff a^^ (cf. Fig. |l|) it is of great interest 
to use the 1/to expansion, implemented as an 
effective theory with the (zero-velocity) HQET 
LagrangianQ (F+V^h = , i^hP+ = V^h 7 P± = 

2 > 

Ckqet = a'^' + ^£i'^+0{^), (25) 
Z:^*^' = MDo + mjtPi,. (26) 
The higher order terms such as 

4'^ - ?h(-iD2-Ba)V'h, (27) 

are to be treated as perturbations to the static 
theory defined by L^'^*'. Then the expansion is 
expected to be renormalizable order by order in 
1 /m (similar to chiral perturbation theory) . Nev- 
ertheless, £hqet contains operators of different 
dimension which mix under renormalization with 
power divergent coefficients ^ a^". Estimating 
them at a given order k in the coupling constant 
expansion, one is left with error terms of the form 

and the continuum limit does not exist. 

This is relevant already at the level of the static 
theory: a linearly divergent (71 = 1) additive mass 
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renormalization is needed. Unlike the case of rel- 
ativistic fermions, there is no (chiral) symmetry 
to fix this term. Thus, a continuum hmit of the b- 
quark mass in the static approximation requires 
a strategy for NP renormalization including the 
power divergent subtraction. A viable approach 
was introduced at last year's conference (5|]. The 
simple idea is again based on the use of finite vol- 
ume. It is illustrated in Fig. 0, which one should 
compare to Fig. |l|: if we choose the box size L 
appropriately, we can both 

• treat the b-quark as a relativistic fermion, 
i.e. simulate QCD with ami, 1 and 

• apply HQET quantitatively which requires 
Lmb » 1.0 

Identifying some suitable observables in the 
Mb 

I — I — I — I — I — I — I — z 

4 8 12 Energy 

[GeV] 



Nf = 4 QCD + HQET 



finite volume latt. 



Figure 7. Matching HQET and QCD for L ^ 0.2fm 

two theories allows to connect the parameters of 
HQET to those of QCD. This solves the match- 
ing problem entirely non-perturbatively (in this 
case the intermediate step through some RGI's 
(Fig. |l|) is not necessary). 

The idea can be implemented in such a way 
that only quantities which are defined separately 
in HQET and in QCD and which are independent 
of the regularization have to be computed. To 
understand how this works, it is best to consider 
the following example. 

3.3.1. The b-quark mass in the static ap- 
proximation 

The renormalization of the static theory re- 
quires (apart from the usual one in the light sec- 
tor) a (linearly divergent) additive mass renor- 
malization. In other words, when we insert m — 



^ These conditions arc satisfied for L 0.2 fm, where in 
comparison to a large volume simulation (L » 2 fm) an 
order of magnitude is gained in the lattice spacing for the 
same L/a. 



mbarc(.9o) into the Lagrangian £fj *, eq. (g6|), 
with suitably chosen TObare(5o), then all energies 
are properly renormalized and have a continuum 
limit. Due to the simple form of the action, the 
dependence of the static propagator on mbarc(5o) 
is explicit and corresponds to a common shift of 
all energies. It is therefore notationally simpler 
and common practice to define all observables 
with m = in eq. (^) and keep in mind that 
renormalized energies are 



■stat,R, 



-Estat + mbarc , with ni = 0. 



(28) 



We now want to use the B- meson mass, me, on 
the left hand side and relate mbaic to the quark 
mass in the relativistic theory, whose renormal- 
ization is under control |ll| to obtain the relation 
between tob and the RGI b-quark mass. 



experiment 



Lattice with aniq <^ 1 



ma =5.4 GeV T{Lo,M) 

rstat(-t'2) — ■ , rstat(il) " — , rstat(-^^o) 

Figure 8. Connecting experimental observables to 
properly renormalized HQET. 



Starting from a static-light correlation function 
in the SF (see for the precise definition), 

/r*(xo) c j2 (^r (^) Ch(y)75Ci(z)) (29) 

y >z 



A Stat 
^0 



Cl 

xq — xq = L 

wc first define an i-dependent bare energy 



(30) 



rstat(i) = -^{do + dt)\n[f]^^\xo)] 
at Xq — L/2 . 

The renormalized energy is 

rstat,R(-^) = rstat(i) + ™baro • 



(31) 
(32) 



10 



The same quantity computed in QCD with finite 
mass of the heavy quark is denoted by T{L, M), 
where the argument M is conveniently taken as 
the scheme independent RGI mass [ID of the 
heavy quark. Matching the two theories at lowest 
order in 1 /M is achieved by the condition 



rstat,R(i) = r(i,Afb). 



(33) 



To the same accuracy, the energy of a B-meson, 
^'stat.R, equals the physical (experimental) mass, 
niB- Thus we have 



4.5 



4 - 



3.5 




Lo r„,(Lo,M) 



_J_ 



5.5 



6.5 M L„ 



Figure 9. Mass-dependence of LoT{Lo, M), in the 
continuum limit (A^f = 0). 



TUB 



= £;,tat + r(L,A/b)-r,tat(i) 

= [Sstat -rst,at(i„)] (34) 
+ [rstat(in) — Tstat (Lo)]+r(Lo,Mb), 



for arbitrary L„. Here the matching eq. ( |33| ) is 
performed at L ~ Lq, which has to satisfy 



LoMb > 1 



(35) 



The two terms in "[ ]" -parenthesis are energy 
differences in the static theory which do not need 
renormalization. Their continuum limit can be 
taken. The entire quark mass dependence is con- 
tained in the last term T{Lo, Mb), defined in QCD 
with a relativistic b-quark. For a range of M, 
the function T{L(), M) can again be computed in 
the continuum limit. The result of a quenched 
study [ pT| is shown in Fig. ^. 

Finally, all ingredients are present to solve 
eq. (HJ) for Afb- The various steps of the pro- 
cedure are summarized in Fig. M. 



One point remains to be discussed. Of course 
the computation of i?stat needs a lattice large 
enough to accommodate a B-meson, say of size 
L > 1.5 fm. In this step one can then not af- 
ford lattice spacings much smaller than 0.05 fm. 
At the same lattice spacing (the same bare pa- 
rameters of the theory), one needs to compute 
rstat(iri) to form the difference £^stat -^stat(^'r^)• 
Therefore also Ln may not be too small. In prac- 
tice, Ln = 2"Lo, Lq « 0.2 fm with n = 2 is 
sufhcient[|l). 

To connect Tsta.tiL2) to Tstat{Lo) it is further- 
more convenient to use 

[r.tat(i2) -rstat(io)]io = (36) 
50-111(^0) -I- i(Tin(wi) with Ui = 5^(ii), 



(the Schrodinger functional coupling) and 

o'miu) = lim YiT^{u,a/L) (37) 

a/L^O 

Ein (u,a/L) = 2L [rstat(2i) - rstat(i)]u=ff2(L) . 



Unnecessary additional scales are avoided by set- 
ting the light quark mass to zero everywhere, ex- 
cept for i?stat ■ There it is set to the strange quark 
mass [ p3| . Correspondingly the experimental spin 
averaged mass tob = m-B^ — 5405 MeV enters. 



0.5 



0.4 



0.3 




0.3 




0.1 0.2 0.3 0.4 0.02 0.04 

a/ro (a/ro)^ 
Figure 10. Lattice spacing dependence of the 
subtracted B-meson energy with 0(a) improvement 
(r.h.s.) and without (l.h.s.). 

The method has been tested in the quenched 
approximation [|l|. Except for E'stat - rstat(-^2) 
all pieces entering eq. (|3^) have been extrapolated 
to the continuum. We here just mention few de- 
tails of this numerical work. 

• 0(a)-improvement was employed in all steps. 
Although the improvement coefhcient entering 
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the static axial current is known only perturba- 
tively, this did not play a prominent role. Its 
effect on Em is very small. 

• The continuum extrapolation of r(Lo, Af) re- 
quires to fix Lq and z = MLq while varying Lo/a. 
The proper values of go to keep g^{Lo) = 2.45 are 
known from pO|jll| | and the bare subtracted quark 
mass rriq of the heavy quark is obtained from 

z — MLq = Lq Z Zm mq{l + bn^aniq) (38) 

The coefficients Z , Zm , are defined in 
imiS^ and have also been computed non- 
perturbatively Q. Nevertheless 5m and Z are 
needed at relatively small go, outside the range 
considered in They needed to be extrapo- 

lated as a function of go, causing additional un- 
certainties, which can be reduced in future work. 

• Values 5.15 < z < 9 were used. In this range, 
r(Lo, M) turned out to be an almost linear func- 
tion of z as naively expected. It could be inter- 
polated in z easily (see Fig. ^). 

The last piece, i^stat — rstat(i2), has now been 
computed for both the unimproved theory (i?stat 
from [Q) and with improvement |^ for various 
lattice spacings. Still the situation, shown in 
Fig. |l^ is not completely satisfactory. The im- 
proved results show some residual dependence on 
the value of the time-separation where -Bstat is ex- 
tracted; this leads to large overall uncertainties. 
The unimproved results generally are somewhat 
higher than the improved ones, but do have large 
errors for a < 0.1 fm. We consider 



t (L2)]Lo- 0.40(4) 



(39) 



which is shown at a = 0, as an estimate of the 
continuum (quenched) result. Setting the scale 
through To and using Lo/ro ~ 1.436/4 the 
solution of eq. (|3^) reads 

(40) 



Afb = 6.96(8) (10) GeV, 



where the larger part of the error originates from 
the uncertainties of Z , Zm , 6m ■ 

3.3.2. Generalizations, Perspectives 

The above example shows that the matching 
in finite volume can be used in practice and for a 
power-divergent quantity. We now explain briefly 
how it can be applied more generally. 



Again consider multiplicative renormalization, 
e.g. of the static-light axial current. In terms 
of a matrix element, X{L), of the bare operator 
defined in finite volume, the renormalized matrix 
element is 



XniL) = Z{go,^^a)X{L)., 



(41) 



with a renormalization constant Z{go, /ia) defined 
in some independent way (here fi and L are not 
related). The desired large volume renormalized 
matrix element is denoted by 

$R = Z(5o,H*- (42) 

In complete analogy to eq. (^^ we then have[] 

*R - $r"* + 0(l/m) (43) 



^stat 



XH-^L^) X^'^*(Li) XH-^Lo) 



xr\Lo) 



o-x(wi) crx(uo) ^r(-^o, M) . 



The replacement of X^^^Lo) by Xn{Lo,M) is 
the matching in small volume. The step scaling 
functions are (Z(go,iJ-a) cancels out!) 



(7x{u) — lim 



j^stat^2i) 



(44) 



These formulae are another explicit realization of 
the strategy Fig. || and are easily generalized to 
include mixing and power subtractions. 

All of the above may be applied also beyond 
the leading order in 1/m. This opens the exciting 
perspective to obtain HQET predictions includ- 
ing the first order 1 / m-correction and completely 
non-perturbative in the QCD coupling. The 1 / m- 
terms are then treated as insertions in the static 
propagator (see e.g. Q). Of course, the match- 
ing of all quantities has to be done consistently 
to order 1/m. 

A potential obstacle on this road to precision 
predictions for B-physics is the generally rather 
bad signal-to- noise ratio of static- light correlation 
functions in large volume ||5^. NRQCD suffers 
little from this problem. Unfortunately in this 

^In the case of mixing, the step scaling functions become 
matrices and also the first factor in the last line can easily 
be given a proper definition. 



12 



theory the dimension five operator iD^)?/'h 
is part of the leading order Lagrangian rendering 
the theory unrenormaHzable. It appears impos- 
sible to formulate NP-renormalized NRQCD in 
such a way that a continuum limit can be taken. 

4. CONCLUSIONS, OPEN PROBLEMS 

Over the last years NP renormalization has 
been developed and applied at a steady pace. Ma- 
jor recent steps have been 

• the application of the MOM-scheme with 
fermions with exact chiral symmetries, 

• the renormalization of the static-light axial cur- 
rent 

• first steps towards the renormalization of 4- 
fermion operators in the SF 

• the application of the SF renormalization for 
Nf = 2 and 

• the development of a strategy to renormalize 
HQET non-perturbatively. 

The latter, proven to work in an example, should 
be developed further, in particular to include 
1/m-corrections. Important questions are in how 
far it is necessary to implement Symanzik im- 
provement and how one can reduce the noise of 
long distance correlation functions. On the other 
hand, no particularly small lattice spacings are 
needed and thus the approach can be used for 
dynamical fermions. 

The computation of asF (p) and the renormal- 
ization of composite operators in the SF-scheme 
have been shown to be applicable in the presence 
of dynamical fermions without excessive com- 
putational requirements. Unfinished parts are 
mainly due to the connection to the low-energy 
region of the theory where the usual problems 
(are the quarks light enough, the volumes large 
enough?) have to be tackled. 

Of course, concerning e.g. asFi^J■)^ one needs 
iVf = 3 dynamical fermions to make better con- 
tact with reality. A first step in this direction 
has been taken by the JLQCD/CP-PACS collab- 
oration who have computed the important im- 
provement coefficient Csw non-perturbatively for 
iVf = 3, using both the Wilson gauge action and 



the Iwasaki one . 

An important open problem is whether the SF 
may be combined with fermions with exact chi- 
ral symmetry. It should be possible to formulate 
the SF with any kind of regularization, since it 
is (formally) defined in the continuum. However, 
no convincing formulation is known as yet. This 
would be important in order to be able to com- 
pute the last factor Winter (ffo j Mmo), eq. (|8|), also 
in such regularizations without going through yet 
another quantity defined with periodic boundary 
conditions or in infinite volume 
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